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How much does the curvature perturbation change after it leaves the horizon, and when should 
one evaluate the power spectrum? To answer these questions we study single field inflation models 
numerically, and compare the evolution of different curvature perturbations from horizon crossing 
to the end of inflation. In particular we calculate the number of efolds it takes for the curvature 
perturbation at a given wavenumber to settle down to within a given fraction of their value at 
the end of inflation. We flnd that e.g. in chaotic inflation, the amplitude of the comoving and the 
curvature perturbation on uniform density hypersurfaces differ by up to 180 % at horizon crossing 
assuming the same amplitude at the end of inflation, and that it takes approximately 3 efolds for 
the curvature perturbation to be within 1 % of its value at the end of inflation. 

PACS numbers: 98.80.Cq |arXiv:1111.6940l 



I. INTRODUCTION 

Our understanding of the early universe has rapidly improved over the last couple of decades. Observations indicate 
that we live in a universe that agrees remarkably well with the cosmological standard model, a universe that is 
homogeneous and isotropic on the largest scales, described by the Friedmann-Robertson- Walker (FRW) spacetime 
and that includes cold dark matter and a cosmological constant. On smaller scales, anisotropics in the Cosmic 
Microwave Background (CMB) and the Large Scale Structure (LSS) are sourced by quantum fluctuations formed in 
scalar fields during an inflationary epoch at very early times. As the process describing the decay of these scalar 
fields into standard matter is not known, we map the power spectrum of the field fluctuations onto the spectrum 
of a conserved quantity when the size of the fluctuations is similar to the size of the Hubble radius, that is at 
horizon crossing, k = aH, (see e.g. Refs. [H |2j for reviews of inflationary cosmology). The relative size of a given 
field fluctuation compared to the Hubble radius (l/aH) is important. Perturbations arising from fluctuations in the 
inflaton scalar field will become fixed, or are "frozen in" , when that scale first leaves the Hubble radius. They will be 
conserved while that mode is beyond the Hubble radius if mapped to a suitable quantity like a curvature perturbation, 
such as those discussed in the following. 

Observational cosmology is entering an era in which the data from observations of LSS and the CMB are becoming 
much more detailed. Only eight years ago, the WMAP team were quoting cosmological parameters to an accuracy of 
about 10% , Ref. [31 • Now, as data sets are improving both in quality and size, the WMAP seven year observations 
(hereafter WMAP7) can constrain these parameters to within a couple of percent, Ref. [1]. With Planck we expect 
to do even better. Hence it is essential that the quantities we wish to study in the early universe are understood, 
theoretically, to this same level of precision. 

It is well known that the curvature perturbations on both uniform density hypersurfaces, C, and on comoving 
hypersurfaces, TZ, are conserved on large scales where gradient terms can be neglected (i.e. in the limit k ^ 0) for 
adiabatic perturbations, Refs. 0(5]. This result follows from the conservation of energy, Ref. |7|. The standard 
approach used to calculate the power spectrum of perturbations after horizon crossing assumes that the limit fc — ^ 
has been reached, see Ref. [5]. However, immediately after horizon crossing the wavenumber will not yet have become 
sufficiently small for this limit to be accurate and gradient terms will still play a role. In fact, in single field inflation 
there will still be some residual non-adiabatic pressure perturbation, i5Pnad, present Ref. [9]. So, even in the absence 
of other sources of ^Pnad, the curvature perturbation will continue to evolve for some number of efolds before settling 
down to its value at the end of inflation. Although it is well known that this evolution continues for a short time 
after horizon crossing, the exact amount of evolution has not been quantified. Exactly how long the evolution will 
last and how big the errors may be if the curvature perturbation is evaluated too early are issues which are yet to be 
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addressed in the literature. 

Much work has been done calculating the power spectra for the curvature perturbations, ^ and TZ, during inflation 
(see for example Refs. pDHTS] and the reviews Refs. [121 HZ]). Analytic studies have to rely on either the slow roll 
limit, or large scale approximations, to make the calculations viable. In these limits, and without anisotropic stress, 
the two definitions are equal up to a sign difference e.g.Ref. |18j . In this paper we will quantify the evolution of the 
curvature perturbation shortly before, during and after horizon crossing. In Section |IV| we will demonstrate how C 
and TZ differ at these times, how much the instantaneous horizon crossing values differ from the values at the end of 
inflation and how long it takes for the quantities to reach these final values. 

As we do not want to rely on the slow roll approximation or the large scale limit, we solve the Klein-Gordon 
equations numerically, having chosen a particular scalar field potential. The linear scalar field perturbations are 
evolved from an initial Bunch-Davies state well inside the horizon until the inflationary expansion ends. Most of the 
results we present in this paper use single field chaotic inflation, the simplest single field inflation model which is in 
agreement with WMAP7. In order to ensure that our results are representative beyond this simplest model, we also 
study a set of more complicated single field models, U = Uq + ^m?ip'^, U = and U = oLp^l'^ . We compare these 

numerical results with the standard analytic solutions for single field inflation. It is important to recognise that the 
analytic solution is derived using the fc — !■ limit but is evaluated using quantities at horizon crossing. Although this 
has been known in the literature for many years it is not often made clear when this mixing of late time solution 
and horizon crossing values is being used. In this paper we highlight the large magnitude of the inaccuracies which 
would result if a naive calculation of the power spectrum at horizon crossing is performed using horizon crossing values. 

Throughout this paper, although we use quantum initial conditions, we only consider the evolution of the scalar 
field perturbations classically. There are issues concerning how and when the quantum-to-classical transition takes 
place but we will not attempt to address these here. For further discussion on these issues, see for instance, Ref. 

[IHIEQ]. 

The paper is organised as follows. In Section [Tl] we define the relevant variables and present the relevant governing 
equations to determine their evolution. In the following section we detail the numerics we use to solve the equations. 
In Section [IV| we present our results, and we conclude in the final section. 



II. EQUATIONS 

In this section we will derive the equations for the curvature perturbation and (5Pnad- Throughout the calculations 
we expand quantities into background and perturbative contributions, and shall consider scalar perturbations about 
a homogeneous, isotropic, FRW background model. To first order this leads to the line element 

ds^ = -c? (1 -f- 20) drf- + 2a^B^,dr]dx' + [(1 - 2^) % -f- 2E^,^\ dx'dx^ , (2.1) 

where a = ci{r]) is the scale factor, dij is the flat background metric, 4> the lapse function, and ip the curvature 
perturbation, Bi and Ei are scalar perturbations describing the shear (cTg = —B + E'), and r] is conformal time, 
related to coordinate time t by adi] = dt (for details on perturbation theory see e.g. Ref. [18], the notation of which 
we follow where possible). Derivatives with respect to conformal time are denoted by a dash. Greek indices, /i, A, 
run from 0, ... 3, while lower case Latin indices, i, j, k, run from 1, ... 3. We shall employ the flat slicing and threading 
(or flat gauge) throughout this paper, where not stated otherwise. 



A. Governing equations 

In this paper we only consider single held inflation, which is governed by the Klein-Gordon equation, Ref. [2], in 
the background, 

^'^ + + a^U^^ = , (2.2) 

where H is related to the Hubble parameter hy H = aH, and prime denotes a derivative with respect to conformal 
time ?7, and the Friedmann equation, 



31/2^ V2 



<Po' + a'C/o) • (2.3) 
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By perturbing the field and using the metric defined above in Eq. (2.1 ), we arrive at the perturbed evolution equation 
in Fourier space, 

5^" + + + 0? + (^2^^,C/,^ + ip'^^ ,^^U^ | = , (2.4) 

where we have chosen the flat gauge. Note that in flat gauge the field fluctuation is also known as the Sasaki-Mukhanov 
variable, Ref. [211 [22]. The Fourier component of the field fluctuation, 5ip{k^), is related to the fluctuation in real 

space, 5ip{x^), by 

= (^/ d^k5v{k^)e^^^' . (2.5) 

In order to evaluate the curvature perturbations we will need to relate the pressure, Po and 5P, the energy density 
Po and 5pi to the scalar held, (p^ and (5<^i, by 

^0 - ^^o' - C/o , Po = ^V'o' + t/o , (2.6) 



and 



5p = ^^^W + f C/,^ - ^ 5^. SP= - (u^ + ^ ^/'i ) 6^ , (2.7) 



see e.g. Ref. for details. 



B. Conserved quantities and 5Pnad 

As pointed out above we map the initial scalar fleld fluctuations, which themselves would evolve after horizon exit, 
onto conserved quantities, that remain constant in the limit fc — ^ for adiabatic perturbations. We focus here on the 
curvature perturbation on uniform density hypersurfaces (, Ref. [5], and the comoving curvature perturbation TZ, Ref. 
[6] . The spectrum of the curvature perturbations can then be used to set the initial conditions for standard Boltzmann 
codes see e.g. Ref. [23j, that are used to calculate the CMB anisotropics. The non-adiabatic pressure perturbation 
<^Pnad is in general not directly observable, but is a source term for the evolution of the curvature perturbations, see 
Ref. [7]. Now we use the equations in the previous section above to flnd expressions for the curvature perturbations 
and the non-adiabatic pressure perturbation (5Pnad- 

The total pressure perturbation is split into an adiabatic and non-adiabatic part as Ref. [24] 
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tSp, 



(2.8) 



where Cg is the adiabatic speed of sound and is defined as Cg = Pq'/Po'- This gives us an expressions for SPnad in 
terms of the scalar field quantities, see Ref. [5S] 
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(2.9) 



The curvature perturbation on uniform density hypersurfaces is defined as 



Po 



(2.10) 



which simplifies if we evaluate the right hand side in fiat gauge to —( = ^^pdat , where for ease of use we drop the 
subscript flat in the following. The comoving curvature perturbation, that is the curvature perturbation evaluated on 
comoving or uniform field slices, is defined by 



(2.11) 
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which simphfies again if the RHS is evaluated in flat gauge to 

curvature perturbations, defined in different gauges, are related by the constraint equation. 



■^S(p. Note, that these two gauge-invariant 



(2.12) 



where = ip + 'Hug is the curvature perturbation in longitudinal gauge, that is on uniform shear hypersurfaces. As 
can be seen from Eq. ( 2.12[ ), C + TZ will become small on super-horizon scales. In the work that follows as well as 
studying how quickly TZ and Q settle to their final conserved values we will also be looking at how quickly after horizon 
crossing these two descriptions of the curvature perturbation agree. 

Combining our results above we finally get for an arbitrary wavenumber k for a single scalar field and without 
imposing slow-roll in terms of the scalar field variables the non-adiabatic pressure 
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2Ua^ 



and similarly the curvature perturbation on uniform density hypersurfaces, 

1 



^ o / 2 
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<^0 



2[/a2 



5(p 



(2.13) 



(2.14) 



and the comoving curvature perturbation, defined in Eq. (2.11 ), was already given in terms of the field (/s. 



C. Analytic Solutions 

Using numerical techniques we can evaluate the expressions above at any time, analytically this is not possible. To 
deduce an analytic expression for the curvature perturbation, such as that used in the popular bN formalism, ^ we 
start from the definition [5] 



2^ 



(2.16) 



By obtaining an exact solution for power-law infiation and taking the fc — > limit, and then making an expansion 
about this solution, we arrive at the following, see Ref. [16] 



Vn{kUt2* - [1 - (2C + l)esR + Ct^sk? 



k=aH 



(2.17) 



where C = —2 -t- ln2 + 7 and 7 is the Euler constant and csr and "qsR are slow roll parameters defined by, 
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VSR 
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(2.18) 
(2.19) 



The expression given in Eq. (2.17) is only valid in the large scale limit, or equivalently a long time after horizon exit. 
However, it must be evaluated exactly when the corresponding mode crosses the horizon. In order to arrive at an 
expression which is valid at all times, we would need to consider the full analytic solution, including log corrections. 



see for instance Ref. [30j. We shall return to this subject in Section IV D 



^ The 5N formalism uses the direct relation of the curvature perturbation to the perturbed number of efolds on large scales, see, for 
instance, Ref. I26I429| . We can relate the linearly perturbed number of expansions to the density by 



dN ^ 1 d'^N ^ 2 



(2.15) 
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The above result holds, to the lowest order in slow roll, for any general potential and if the slow roll parameters 
are assumed to be very small this simplifies to 

-PnikUi*^ ff,, ■ (2.20) 

(27r)^(^Q k=aH 

For further discussions on the analytic treatment of curvature perturbations close to Horizon crossing, see for example 
Ref. [311 122] 

III. NUMERICAL SETUP 

To set up the numerical system we have followed the work done by Salopek et al. Ref. [33j . We select a finite range 
of k modes which will cover all the modes which have been observed in the CMB. The WMAP team have released 
results corresponding to the range k G [3.5x10"'*, 0.12]Mpc~^ so we will consider a similar range below. We use the 
number of efolds, M ~ log{a/ainit), as our time variable instead of conformal time, where Cinit is the value of a at 
the start of inflation and is evaluated by setting a — I today and using the background run, assuming instantaneous 
reheating. 

The initial conditions for the background system are selected depending on the choice of potential. For the potential 
with m = 6.32 x IO'^Mpl, we use initial conditions ipo = 18MpL and ipo,/^ — — O.lMpL, see Ref. |34j . 
Following Ref. |33j we set the initial conditions for each k mode a few efolds before horizon crossing when the initial 
time A/init {k) is such that 

aH\init ^ ^ 

At early times we assume the Bunch-Davies vacuum and we get the following initial conditions for modes well inside 
the horizon, see Ref. [33H35] : 

S^mU = -— ^e-^^''fl + *A') , (3.2) 
aMpLV2fc V aHJ 

where 77 = —{aH{l — €h))~^ and ch — ——h- ■ 

The behaviour of all the modes over the scales we consider are similar, so for clarity only fci, k2 and k^ modes, 
given below, have been shown on the graphs. It is worth noting that k2 is the WMAP pivot scale. 

fci = 2.77 X IQ-^^Mpc-^ = 7.28 x IO^^^^pl (3.3) 
fca = 2.00 X 10"^Mpc"* = 5.25 x 10~®"Mpl {WMAP) (3.4) 
fcg = 1.45 X IQ'^Mpcr'^ = 3.80 x 10"^*7\fpL (3.5) 

The results presented below are for the potential U — ^m^ip^ with m — 6.32 x lG~^A/pL, this value has been 
chosen such that V-nik) = 2.45 x 10^^ at the end of inflation for the WMAP pivot scale. We also obtained results 
for three additional potentials, U = Ua + ^rn^tp^, U = jXf'^ and U — aip'^^^, using parameters and initial conditions 
specified in Ref. [331 ES] ■ All the numerical results have been verified with a second numerical program to ensure 
their accuracy, see Ref. [36] . 

IV. RESULTS 

Having setup the governing equations and the numerical formalism in the previous sections, we can now turn to 
answering the questions we raised in the introduction. Here we present our results which quantify the difference 



between ^ and TZ, the magnitude of error incurred if we would use Eq. (2.161 evaluated at horizon crossing instead of 



Eq. (2.171, and the length of time taken for the power spectra to settle to their values at the end of inflation. 

Before discussing our results in detail it is worth considering the evolution of the power spectra in general. The 
evolution, or indeed conservation, of the curvature perturbations has been studied in detail in the past, see Ref. 
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[5HZ1 inZ] and, as expected, we find that some time shortly after horizon crossing there is no longer any appreciable 
evolution in either the power spectrum of C, P(j(fc) or the power spectrum of TZ, V-jzik), see Fig. 1(a) The values of 
these power spectra converge very quickly onto the same conserved value. This behaviour is also supported by the 
second graph in Fig, 



1(b) which shows the non-adiabatic pressure perturbation, SPnad- 
that during and after horizon crossing the size of JPnad drops sharply towards zero, 
perturbation is directly related to the curvature perturbation (on large scales), see Ref. 

C' OC (5P„ad , 



We find, again as expected. 
The non-adiabatic pressure 

(4.1) 



and hence the rapid decrease in the non-adiabatic pressure perturbation causes the curvature perturbations to settle 
onto a conserved value. 

However, it can also be seen from Fig. 1(a) that there is some evolution of the power spectra immediately after 
horizon crossing. This is well known in the literature where the phrase 'soon after horizon crossing' is commonly used 
to refer to the time at which the power spectra have settled down. In the sections that follow we will be looking at 
this evolution in more detail and in particular quantifying exactly how soon after horizon crossing the power spectra 
reach the final value and how different this is to the horizon crossing values. 



0.001 
0.0001 
^ le-05 
&H le-06 
le-07 
le-08 
le-09. 



H.C. 



10 20 30 40 50 60 70 80 

(a) The evolution of the power spectrum of 'P(;{k) and TZ, 

Vizi^) plotted against the number of efolds, M ■ Both 
power spectra stop evolving shortly after horizon crossing, 

however a short period of evolution immediately after 
horizon crossing is visible as is a difference between (fc) 
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(b) The evolution of (5Pnad is plotted against the number of 
efolds, M ■ i5Pnad rapidly drops towards zero after horizon 
crossing. We apply a cut off to the graph at 10"^'', beneath 
which numerical noise dominates. 



FIG. 1: Both the graphs above are plotted for three different k modes. The black line (left) is fci, the blue line (middle) is ^2 
the WMAP pivot scale, and the red line (right) is k^. 



A. How do C and 7^ differ? 



The curvature perturbation on uniform density hypersurfaces, C, and the comoving curvature perturbation, TZ are 
often used interchangeably. Although they have different definit ions, see Eq. (2.10) and Eq. (2.11), it is well known 
that on large scales they are equivalent, as can be seen from Eq. (2.12). This equivalence is however only strictly true 

In Fig. 



la 



we see that there is a difference 



in the large scale limit, and on smaller, finite scales this is not the case 
between the two curvature perturbations near to horizon crossing. Three different fc-modes are plotted throughout 
there evolution, from deep within the horizon through horizon crossing (indicated by the dotted lines), until the end 
of inflation. 

Fig. [2] shows that C, is as much as 20% larger than TZ at horizon crossing and remains significantly larger for at 
least a couple of efolds. This highlights the importance of making explicit the choice of curvature perturbation when 
carrying out calculations close to horizon crossing. We also note that C and TZ take slightly different amounts of time 
to settle down after horizon exit, as detailed in Section [IVC[ 
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FIG. 2; The ratio of ^ and TZ is plotted against the number of efolds, N ■ Evolution is visible for a short period after horizon 
crossing, during which time C, is noticeably larger than 7?.. (Black line, left: fci. Blue line, middle: k2, Red line, right: fcs) 



B. What is the error incurred by using horizon crossing quantities? 




FIG. 3: The percentage difference between Vik) evaluated, numerically, at the end of inflation and the V{k) obtained at each 
time step. There is as much as 180% difference in the values of 'P(^{k) and 100% difference in the values of Vizik). Both the 
graphs above are plotted for three different k modes. The black line (left) is k\, the blue line (middle) is ^2, the WMAP pivot 
scale, and the red line (right) is fea. 



As mentioned above it is well known that the values of the curvature perturbation power spectra are not the same 
at horizon crossing as they are at the end of inflation. However, it is not clear exactly how much of an error would 
be incurred if one would use the horizon crossing values instead of the correct values at the end of inflation. In Fig. 
[3] we can see exactly how different the power spectra are at horizon crossing compared to the values they take at 
the end of inflation. In Fig. 3(a) Vn{k) is 100% larger at horizon crossing and in Fig. 3(b) 'P(^{k) is 180% larger. 
In fact if we use an analytic approximation that does take into account the pre horizon behaviour we can see the 
100% difference in the power spectra. Near horizon crossing the scalar field's wavefunction (-0) is approximately 
proportional to (1 — ikrj)e^^'^ . At horizon crossing \krj\ — 1 and l-i/jp (x |1 — ikrj\ = 2. A few e-folds later, \kr]\ « and 



cx 



|1| = 1. This is a drop of 1. Hence the power spectra at horizon crossing is expected to be roughly 100% larger 



than that at late times. The factor of two difference in Vnik), has previously been found analytically, for example 
Ref. |19| This is, however, a large difference and might impact on calculations if V{k)k=aH is used to approximate 
the value of the power spectra at the end of inflation. 
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C. How quickly does the power spectra reach its final value at the end of inflation? 




FIG. 4: These graphs show how many efolds one needs to wait for 'P(fc) to be within 10% and 1% of 'P(fc) at the end of inflation 
for a given value of k (measured in Mpc~^). Note that this is independent of k. Table |l] gives the exact number of efolds for 
various potentials. 



The expressions "soon after horizon crossing" and "a few efolds after horizon crossing" are often used in the 
literature to determine when it is reasonable to evaluate the power spectra, so that one can be confident that they 
have the same value as at the end of inflation. Now that we have established that the values at horizon crossing can 
be as much as 180% greater than those at the end of inflation we can study exactly how long after horizon crossing 
we must wait until the power spectra converge onto these final values. Fig. [4] shows how many efolds one needs to 
wait for the power spectra to be within 10% and 1% of the value they have at the end of inflation. It is clear from 
these graphs that this is independent of k and is in the region of "a few efolds" . Table |T] shows the number of efolds 
it takes to be within a flxed percentage of the flnal value, and also shows data for the three additional potentials 
we investigated. The choice of potential makes very little difference to the result. For every potential considered it 
took 1.40 - 1.44 efolds for V(;{k) to be within 10% and 2.56 - 2.86 to be within 1%. It took 1.14 - 1.20 efolds for 
Vizik) to be within 10% and 2.30 - 3.21 efolds to be within 1%. As we highlighted in section |l] observational data 
will soon be constraining observables to within a percent, so this is at least the accuracy we would like to be able 
to evaluate quantities to. Our results show that to be within 1% of the correct power spectra value at the end of 
inflation, evaluating VTi{k) approximately 3.2 efolds after horizon crossing and evaluating V^^ 2.9 efolds after horizon 
crossing, would be sufficient, for all the potentials we studied. 
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TABLE I: The values in this table represent how many efolds after horizon crossing it takes for the power spectrum to be 
within a fixed percentage of the power spectrum at the end of inflation. 



D. Highlighting the magnitude of the error made in using the incorrect analytic expression 



In most analytic calculations an expression for the power spectrum is derived using the /c — > limit but is evaluated 
using quantities at horizon crossing, as given in Eq. (2.17). As we mentioned in sectional it is often not made explicitly 
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(a)The lines in this graph are shown for the mode k2 which 
is the WMAP pivot scale. 
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(b)The lines in this graph are shown for the mode k2 which 
is the WMAP pivot scale (left hand lines) and fca (right hand 
lines) . 



FIG. 5: The evolution of the two types of curvature perturbation, Vizik) and V^ik) are plotted against the number of efolds, 
A/". These numerical solutions are compared to the correct analytic solution with and without slow roll corrections (Eq. (2.171 
and Eq. ( |2.20[ ) and the naive approximation of evaluating the power spectra at horizon crossing. 



clear in the literature when the mixing of late time solution and horizon crossing values is being used, even when this 
is well understood by the authors. In Fig. [5] we compare both this correct analytic solution and the naive calculation 
of the power spectrum at horizon crossing performed using horizon crossing values with the numerical solutions. In 
Fig. [Sj as expected, the correct analytic solution gives a very good estimate to the full numerical solution. Even the 
analytic solution without the slow roll correction, given in Eq. (2.20), is a very good estimate to the full solution, in 



fact the error in not including the slow roll corrections for the WMAP scale is only a slight underestimate of 0.38%. 
However, when we compare the numerical solution to Vnik) evaluated at horizon crossing we find, as shown in Fig. 
that there is a 100% error in our answer. 



3(a) 



3(b) if we evaluated the VQik) 



As we have shown earlier in Fig. 
at horizon crossing the error would be even higher at 180%, this difference in the two power spectra can be seen 
clearly in Fig. [Sj This again highlights the importance of establishing that the correct analytic expression is used in 
calculations, and when different calculations and results are compared. 

Another possible source of error would be to use the analytic expression given in Eq. (2.17) above, but not to 



evaluate it at horizon crossing. If one were to evaluate this expression 'some efolds after horizon crossing', one would 
underestimate the amplitude of the power spectrum. This corresponds to following the red line in Fig. [5(b)| For 
example, evaluating Eq. (2.17) 4 efolds after horizon crossing would incur an error of 15%. Furthermore, evaluating 
the power spectrum at later and later times will increase the error. Lastly, it is worth noting that although the 
analytic and numerical expressions agree in the large scale limit, they do not agree with each other shortly after 
horizon crossing. As shown in the section above one must wait at least 3.2 efolds for these two values to agree. This is 
particularly important if there is a second phase of evolution caused for example by a second scalar field which starts 
to dominate during these three efolds, see for instance Ref. [5T] . 



E. The Spectral Index 



Using observations we can gain information about the power spectrum of curvature perturbations, which in turn 
allows us to constrain cosmological theories. Many observational results are given in terms of a few observables which 
can then be compared directly with predictions given by theories, one such observable is the spectral index. The 
spectral index describes the scale dependence of the power spectrum of the curvature perturbation and is defined by 

m 



dln7\(k) 
dink 



^ dlnVnjk) 
~ dink 



(4.2) 



As an example of how the results presented in this paper will impact on particular observables we consider the 
spectral index in more detail, we present results here using chaotic inflation. Figure [6] shows, as expected, that just 
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FIG. 6: The evolution of the spectral index for the two curvature perturbations, n-jz and n^, are plotted against the number of 
efolds, A/". The lines in this graph are shown for the mode k2 which is the WMAP pivot scale. 



like the curvature perturbation the spectral index continues to evolve for a few efolds after the mode has crossed 
outside the horizon. Evaluating the spectral index naively at horizon crossing gives a results of n-ji « 2 and rt^ ~ 2.3, 
which is an error in both cases of more than 100%. We find that in order for and n-ji to be within 1% of their 
values at the end of inflation, we should evaluate them at least 2.91 and 2.66 efolds after horizon crossing, respectively. 
These values are similar but slightly less than the number of efolds it takes for the power spectrum to be within 1% 
of it's final value, see Table 1. 



V. CONCLUSIONS 



In this paper we have quantified the evolution of the curvature perturbations after inflation and highlighted possible 
errors which can occur. As we are entering an era where we can hope to constrain cosmological parameters to within 
a percent using the observational data from e.g. Planck, it is of particular importance that these errors are both 
minimised and quantified. 

After presenting the evolution equations for the scalar field, we gave the expressions for the curvature perturbations 
we consider. We then outlined the numerical methods used to solve these equations and gave details of the various 
models and the initial conditions used. In presenting our results, we found that despite C and TZ being equivalent very 
far outside the horizon, the difference between \Tl\ and |C| at horizon crossing can be as much as 20%. We also found 
the error in evaluating the power spectra numerically at horizon crossing instead of either using the correct analytic 
expression or the full numerical solution at late times can be as much as 180% for P({k) and 100% for Vn- Lastly we 
showed that if one wanted to evaluate the power spectra without the use of the analytic expression Eq. (2.17), one 
would need to wait at least 3.2 efolds to ensure the answer for V-nik) is correct to within 1% of the value at the end 
of inflation, and one would need to wait at least 2.9 efolds to ensure the answer for is correct to within 1%. There 
was no significant difference to these results when we considered the three additional single field models presented at 
the end of Section HTll 

In this paper we point out that there is a difference between analytic and numerical expressions close to the horizon. 
The numerical results are the instantaneous values of the power spectrum and spectral index at horizon crossing, not 
the expected late time values. These instantaneous results, while not of observational significance, are useful in 
many ways, including as initial conditions for other analytical or numerical schemes which operate purely outside 
the horizon. Firstly, if we are interested in the late time values we should not take the 'naive' numerical approach 
of evaluating these at horizon crossing. Unlike the case when using the analytic expressions these results will not be 
close to the correct answer. Secondly, if we are interested in the instantaneous values at or close to horizon crossing, 
for instance when developing codes which rely on this information, the normal analytic expressions will not give the 
correct answers as they are no longer valid and you need to use numerical methods. 

In this paper we have only studied single field inflation models, in which the non-adiabatic pressure decays rapidly. 



see Fig. 1(b) This is no longer the case in multi- field systems, as a recent work has detailed, Ref. [51]. It will 
be interesting to repeat our analysis for more complicated models where superhorizon evolution of the curvature 
perturbation is expected. This would include, for instance multi-field inflation or "ultra slow roll inflation" , see for 
instance Ref. [131 [SI]- Since in this paper we focus on the single field case, we postpone a study of these cases for 
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future work. 

We have chosen the end of inflation as a natural end point of our calculations. After the end of inflation the 
inflaton is assumed to decay into the standard matter fields during reheating, the detailed mechanism of which is as 
yet unclear. Also, after reheating we have a multi- fluid system, and hence automatically SPnad 7^ 0, see Ref. 39J, 
which means that the curvature perturbations are no longer conserved (on any scale) . Consequently the values at the 
end of reheating might no longer be the same as those at the end of inflation. 

In conclusion, we have highlighted the fact that the different curvature perturbations do evolve differently imme- 
diately after horizon exit. Confusion of the different curvature perturbations can introduce additional errors when 
comparing theoretical results with observations, which can easily be avoided. 
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